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THERMAL  CALCULATION  OF  A FREEZING  TUBE  IN  WATER-PERMEABLE  GROUND 


Doctor  (Technical  Science)  B.  V.  Proskuryakov 


1 . Foreword 

The  problem  of  ground  freezing  has  a theoretical  and  an 
experimental  history.  The  problem  of  freezing  homogeneous  ground  by  a 
single  tube  or  by  a series  of  tubes  has  been  fairly  well  clarified  by  both 
theory  and  experience.  This  cannot,  however,  be  said  about  freezing  ground 
in  the  presence  of  filtration  water.  With  the  exception  of  several  attempts 
to  solve  this  problem  (see,  for  instance  the  article  by  docent  Chamyy 
"Soviet  Metropolitan1',  No.  4,  1940)  which  introduced  a number  of  assumptions 
and  roughly  outlined  the  process,  we  simply  find  no  rigorous  solution  to 
this  problem,  and  are  compelled  to  note  that  the  literature  lacks  even  data 
concerning  experimental  studies  of  this  problem.  However,  in  a number  of 
cases  practitioners  encounter  the  necessity  to  freeze  sections  of  ground 
which  are  sometimes  highly  water-permeable. 

In  this  case  it  is  impossible  to  evaluate  by  calculation  neither  the 
required  power  of  the  freezing  devices,  cor  to  select  an  economical  interval 
between  the  tubes,  nor  to  establish  the  freezing  times,  and  the  problem 
has  to  be  solved  during  production  by  trial  and  error,  causing  superfluous 
expenditures  to  be  made  on  freezing  and  losing  time  in  experimenting  during 
the  process  of  carrying  out  the  operations.  It  must,  however,  be  noted  here 
that  there  is  practically  no  experience  in  freezing  ground  under  these 
conditions.  The  goal  of  our  study  is  to  establish  means  of  solving  this 
problem. 

2.  Calculating  the  Freezing  of  Water-Permeable  Ground  by  Means  of  a Single 
Tube  (Steady-State  Conditions) 

Let  us  assume  that  a freezing  tube  of  radius  r has  a surface  temperature 
T which  is  kept  constant  throughout  the  freezing  period. 

This  tube  is  surrounded  by  the  filtration  flow,  the  direction  of  which 
is  perpendicular  to  the  tube  axis  and  the  velocity  of  which  when  learning  the  tube 
is  equal  to  v.  The  temperature  of  the  ground  and  consequently  of  the  filtration 
flow  when  approaching  the  tube  is  equal  to  T^. 

The  ground  will  continue  to  freeze  until  the  quantity  of  heat  emitted 
from  the  surface  of  the  freezing  column  of  ground  to  the  tube  is  equal  to 
the  quantity  of  heat  transmitted  to  this  surface  from  the  filtration  flow. 

When  these  quantities  of  heat  are  equal,  freezing  will  stop  and  a state  of 
equilibrium  will  be  established  in  which  the  power  of  the  freezing  mechanisms 
will  be  expended  only  to  cool  the  filtration  flow  if  the  dimensions  of  the 
column  of  frozen  ground  are  kept  constant.  We  will  examine  in  this  way  the 
problem  under  the  conditions  of  established  heat  conditions  specifically  for 
the  case  where 

Q2  = Qr 
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where  Q_,  - is  the  quantity  of  heat  extracted  from  the  filtration  flow  at 
the  boundary  between  the  thawed  and  frozen  ground; 

- is  the  quantity  of  heat  taken  from  the  surface  of  the  frozen 
ground  (zero  isotherm)  from  the  freezing  tube. 

In  this  case  the  problem  is  stated  as  follows: 

"At  a given  velocity  v of  a filtration  flow  which  has  temperature  T^. , 
find  the  maximum  radius  Rq  of  a frozen  column  of  ground  if  the  radius  1 

of  the  freezing  tube  r^,  and  the  temperature  on  its  surface  T are  known." 

In  solving  this  problem,  we  are  first  obligued  to  state  that  the 
propagation  of  heat  within  the  boundaries  of  the  frozen  portion  of  ground 
and  within  the  limits  of  the  thawed  ground  penetrated  by  the  filtration  flow 
obey  different  laws.  In  the  first  case  heat  propagates  only  due  to  the 
thermal  conductivity  of  the  ground,  and  in  the  second  case  heat  is  transferred 
from  point  to  point  by  the  filtration  flow  as  well. 

Thus,  we  are  obligued  to  examine  two  problems:  One  concerns  the 

propagation  of  heat  in  the  frozen  portion  of  ground  and  the  second  concerns 
the  propagation  of  heat  in  the  thawed  ground  in  the  presence  of  a filtration 
flow,  or  as  we  will  refer  to  it  for  the  sake  of  brevity,  the  internal  problem 
and  the  external  problem. 

In  solving  the  problem  we  will  consider  that  the  frozen  column  of 
ground  is  a circular  cylinder  which  is  eccentrically  located  relative  to 
the  freezing  tube.  This  assumption,  as  experience  shows,  can  be  accepted 
without  introducing  any  large  errors. 

A.  The  Internal  Problem 

Problem  statement.  We  will  designate  the  following: 

1.  The  radius  of  the  frozen  column  of  ground  or,  what  is  the  same 
thing,  the  radius  of  the  zero  isotherm  Rq. 

2.  The  outer  radius  of  the  freezing  tube  r^. 

3.  The  temperature  on  the  surface  of  the  freezing  tube  T. 


4.  Eccentricity  e of  the  zero  isotherm  with  radius  R 

o o 

and  oP  the  isotherm  with  radius  r^. 

5.  The  heat  conductivity  factor  of  the  frozen  ground 

Find:  1.  The  total  amount  of  heat  extracted  from  the  zero  isotherm. 


2.  The  distribution  of  heat  extracted  from  the  zero  isotherm. 


In  view  of  the  fact  that  the  heat  in  the  frozen  ground  propagates 
only  due  to  thermal  conductivity,  we  may  consider  that  in  this  case  the 
Fourier  equation  is  applicable: 


dt  / d-t  d-t  ’ 
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where 

a - is  the  temperature  conductivity  factor; 
t - is  the  temperature  at  any  point  of  the  frozen  ground; 
t - is  time; 

x,  y - are  the  coordinates  of  a point  with  temperatures  t. 

In  view  of  the  fact  that  we  are  examining  established  thermal 
conditions,  at  which  the  temperatures  do  not  vary  in  time,  we  can  write  the 
Fourier  equation  as  applied  to  our  case  in  the  following  form: 


dxi  °y- 


The  solution  to  this  equation  as  applied  to  the  given  problem 
conditions  is  a function  which  is  known  in  the  theory  of  a two-charge 
electrical  field  (a  dipole  field) : 


|InV  -j  - (y  + «)- — In^x-  — (y-tf/-  — C,  J . 


The  position  of- the  coordinate  axes  and  the  accepted  designations 
are  shown  in  Figure  1. 

Differentiation  readily  demonstrates  that  function  (3)  satisfies 
thermal  conductivity  equation  (2)  and  consequently  is  a partial  integral  of 
it . 

Simple  transformations  convert  function  (3)  into  the  following  form: 


:c:‘z±±LYd.  
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The  form  of  this  equation  makes  it  possible  to  state  that  the 
isothermic  curves  are  a family  of  circles  with  radii: 


r »-r-!  - 7-  -.\ 
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and  with  centers  located  on  axes  y,  while  the  ordinates  of  the  circles 
centers  are  equal  to 


e-C,t  - 2C,  j 

>’c=ti  7?,7-^rr  ■ 


By  solving  the  last  two  equations  together,  we  find 


a-  ■=>'  - r- . 


(6) 


(7) 


We  will  determine  the  values  of  the  integration  constants  from  the 
following  condition:  when  y^  = h,  x = 0 we  should  have  r = r^,,  and  t = T. 

Then  from  (5)  and  (6)  we  have 


CJ~t 


- :n 
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From  the  second  condition  for  the  zero  isotherm,  we  find: 

when  y = H and  x = 0,  we  should  have  r = R and  t = 0, 
c o 
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(9) 


By  substituting  the  value  of  C2  from  (9)  into  (8)  we  obtain  an 
expression  for  the  integration  constant  C^: 


C — - 
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In  view  of  the  fact  that  because  of  the  nature  of  the  problem  we 

usually  do  not  know  the  value  of  h,  a and  H,  it  is  necessary  to  replace 

these  values  with  r_,  R and  e . 

To  o 

Using  the  conditions  of  the  problem  and  Figure  1,  we  can  write: 


_ r r r 

P ri / 


(11) 


However,  equality  (7)  makes  it  possible  to  establish  the  following 
functions : 


/ 


Combined  solution  of  equations  (12)  and  (11)  yields: 


> !-r  r,  _ 7/ 
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from  which  we  obtain  the  value  of  a: 


G-’V 

f 


(14) 
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From  (12)  and  (14)  we  find  the  desired  values: 


[(«s-4)-<r , 1 0: 


eo~r(Ro~~rr) 


After  obtaining  the  field  of  isotherm  expressed  by 

e , R and  r^,  we  find  the  value  of  the  heat  flow  directed  from  the 
o o T 

zero  isotherm  to  the  surface  of  our  freezing  tube. 


y = 0. 


line : 


For  this  purpose  we  utilize  equation  (3)  related  to  ordinate 


As  is  evident , in  this  case  the  isotherm  is  a straight 


>-°c7  ‘ U'J 

The  quantity  of  heat  which  passes  through  this  isotherm 
and  consequently  through  any  other,  in  particular  the  zero  isotherm 
(since  the  thermal  conditions  in  our  case  have  been  established)  can  be 
expressed  by  the  following  equation: 


Q ==->-, 


We  find  the  general  expression  for  the  gradient  ^ from 
equation  (3) : 


jy  ,v-  -f  (y  + <i)J 
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Hence , when  y = 0 the  temperature  gradient  will  be  equal  to: 


<).'  ; C,  'l  ) 
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By  substituting  the  value 


from  equation  (20)  into  equation 


(18)  and  by  integrating  the  latter,  we  find: 

-C  { * dx-  ■ 

f C,  (.t-  -rd-)  Ci 


This  equation  can  be  rewritten  as  follows  if  C is  expressed  on 
the  basis  of  (10) : 


II  , / li  , \ 
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The  values,  however,  h and  H can  be  replaced  through  (15)  and 
(16),  after  which  we  find  the  final  value  for  Q: 
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This  solution  is  also  the  solution  to  the  first  portion  of  our 
set  tasks:  Finding  the  total  quantity  of  heat  extracted  from  the  zero 
isotherm.  However,  expression  (23)  contains  a value  as  yet  unknown 


e^  which  has  to  be  determined.  For  this  purpose  we  evaluate  the  distribution 


of  the  heat  extracted  from  the  zero  isotherm.  We  can  state  that 


this  distribution  depends  on  the  eccentricity  eQ  of  the  zero  isotherm 


with  regard  to  the  isotherm  of  T (on  the  surface  of  the 
freezing  tube),  and  the  axes  of  symmetry  is  the  y axis. 


We  will  call  the  heat  losses  through  half  (a  semi-circle)  of  the 
zero  isotherm  directed  towards  negative  values  of  y (or,  when  it  is 
the  same  thing  in  our  case,  towards  the  eccentricity)  Q . We  will  attempt 
to  find  the  relationship 


which  will  also  characterize  the  non-uniformity  of  the  distribution  of 
heat  flow  from  the  zero  isotherm. 


Figure  2. 

Value  is  determined  by  integrating  equation  (21)  within  the  limits 
from  the  flow  line  going  from  point  on  the  zero  isotherm  with 
coordinates  x = -R  , y = H to  a current  line  corresponding  to  point  M.,  on 
the  zero  isotherm  with  coordinates  x = Rq  and  y = H (see  Figure  2). 

Designating  these  integration  limits  as  -X  and  +X,  we  will  express 
them  through  values  already  known  to  us.  Returning  to  Figure  2,  we  see 
that 


(25) 

(26) 


we  find : 


By  integrating  equation  (21)  within  the  interval  -X  to  +X, 


i‘  dx  =~^L( arctan^  —arc  tan— -^|  = ” Cretan— ; 


Q = “ I — =11 
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By  replacing  X in  this  equation  with  its  value  according  to 
equation  (26),  we  obtain: 
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(28) 
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however,  since  a2  = H2  - R , then  we  can  rewrite  equation  (28)  in  the 
following  form: 


'21  f'  - / a\  2X{  j r («3‘ 
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By  substituting  into  equation  (24)  the  values  we  have  found  for  Q and 
Q,  we  find:  e 


E = arc  tan  — 


This  expression  (30)  is  not  by  itself  sufficient  to  find  the  value 
e . It  is  also  necessary  to  examine  the  non-uniformity  of  the  distribution 

of  the  heat  flow  passing  through  the  outline  of  the  zero  isotherm 
based  on  the  conditions  of  heat  emission  by  the  cylinder  to  the  filtration 
flow.  If  in  this  process  it  is  impossible  to  express  the  value  E to  known 
values,  then  the  value  of  e^  will  be  determined  from  equation  (30)  and 

the  problem  will  be  completely  solved. 

Let  us  now  examine  heat  exchange  between  the  cylinder  and  the 
filtration  flow  moving’ around  it. 

8.  External  Problem 

The  solution  of  this  problem  amounts  to  determining  the  heat  losses 
of  a filtration  flow  moving  around  a cylindrical  column  of  frozen  ground. 

As  in  solving  the  first  portion  of  the  problem,  we  will  assume  that 
the  thermal  conditions  have  been  established,  i.e.  , we  will  examine  the  boundary- 
equilibrium  thermal  state  of  the  system. 

In  this  case  the  differential  equation  of  heat  conductivity  may  be 
presented  in  the  following  form:1 


u I — • 


i i I'i 
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•See  F.  Frank  and  R.  Mizes,  Differential  and  Interval  Equations  of 
Mathematical  Physics,  ON’TI  , 1937,  page  675. 


where 


a - is  the  thermal  conductivity  coefficient  of  the  water-permeable 


X 


X is  the  heat  conductivity 
wg 


giound  equal  to  wg  (where  in  turn 
C 

wyw 

factor  of  ground  saturated  with  water  and  Cw  and  yw  are 
respectively  the  heat  capacity  and  volumetric  weight  of  water); 


v and  v - are  velocity  projections  of  the  filtration  flow  rate  onto  the 
x v . . 

coordinate  axes. 

The  method  of  integrating  equation  (31)  for  the  case  where  a 
cylinder  of  infinite  length  is  surrounded  by  a flow  has  been  proposed  by 
Bussineskoe . 1 


We  will  replace  the  variables  x and  y by  c,  and  n with  the  aid  of 
the  following  conversion  formulas: 


di_dr 

oy  dx  ' ’ ux  Oy 


(32) 


where 

n - is  the  velocity  potential  which  is  constant  along  equipotential 
lines ; 

5 - is  the  "current  function",  which  is  constant  along  the  current 
line . 

The  system  of  current  line  and  equipotential  lines  is  orthogonal  and, 
as  is  known  from  hydrodynamics,  satifies  the  equation: 


<H  Or  v : 

0.x  Ox  ' oy  Jy 

Examining  the  left-hand  portion  of  equation  (31)  which  is  in 
parentheses  and  which  utilizes  the  relationship  between  the  variable  E,  and 
n and  the  previous  values  of  x and  y,  (32)  makes  it  possible  to  write  the 
following  as  the  result  of  simple  differentiation: 


O r , (Of  __  j_  o-t 

Ox-  Oy  - x 1 •*'/  \o  j-  0 i,- ' 


(33) 


1 .Analytical  Theory  of  Heat,  Paris,  1905. 


With  allowance  for  (32),  however,  the  right-hand  portion  of  equation  (31) 
acquires  the  following  form: 


ut 

0rt 


(34) 


These  last  two  equations  are  rewritten  with  allowance  for  equation 
(31)  in  the  following  form: 


.$‘V 

01 

i/r-J 

0 r. 

(35) 


If  the  axis  of  the  cylinder  passes  through  the  origin  of  the  coordinates, 
then  the  equations  of  the  current  lines  and  equipotentials  will  be: 


( Pry  \ 

i R ‘ .v 

r.=rul-V\x  + -r-:  , 
\ a- 


(36) 

(37) 


where  V - is  the  velocity  of  flow  in  infinity,  where: 


^ ' ^-inf  ’ ancK'vinf 


(38) 


We  select  the  values  £ and  n such  that  when  /x2  + v2  = R 

oo  '0 

£ becomes  0;  thus  £q  = 0 , and  when  x = Rq  and  y = 0 we  can  approximately 

assume  n = 0,  from  which  n = 2R  V. 

o o 

At  these  values  of  £ and  y equations  (36)  and  (37)  are  rewritten 
in  the  following  form: 


(39) 


(40) 


As  applied  to  the  conditions  of  the  problem  to  be  solved  by  us, 

32t 

it  is  possible  to  ignore  the  value  ^ — r [in  equation  (35)]  in  comparison  with 


assuming  that: 


(41) 
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Then  equation  (35)  acquired  the  form: 


vr 
i)  r 


(42) 


Equation  (42)  should  be  integrated  under  the  following  boundary 
conditions : 

On  the  body  surface  when  £ = 0 t 

when  n = 0 t 

and  everywhere  in  infinity  t 

Whenever  discussing  the  method  used  to  integrate  equation  (42) 
under  boundary  conditions  (43)  and  relying  on  the  special  guides  of 
interest  to  us,1  we  will  cite  only  the  final  expression  for  evaluating 
the  quantity  of  heat  which  the  flow  surrounding  our  frozen  column -cylinder 
loses  (per  one  linear  meter  of  cylinder  length): 


= 0 
= T 

= T 


fi 

fi 


(43) 


n= -l  i -'wg^w 


Tfi.. 


(44) 


where 


is  the  heat  conductivity  factor  of  the  ground  saturated  with 
water-; 


Yw  - is  the  weight  of  a unit  of  water  volume; 


Cw  - is  the  heat  capacity  factor  of  the  water; 


T..^  - is  the  temperature  of  the  filtration  flow  on  the  approach 
to  the  cylinder  (in  infinity). 


In  this  latter  equation  value  corresponds  to  x = Rq  and  y = 0, 
and  therefore  from  (40)  we  obtain: 


(45) 

When  x = -R  and  y = 0,  we  find: 
o 

=4.-V-'-  (46) 


‘For  instance,  Frank  and  Mizes.  Differential  and  Interval  Equations  of 
Mathematical  Physics,  ONTI , ly37,  page  676. 


By  substituting  the  values  of  n and  r\  obtained  from  (45)  and  (46) 
into  equation  (44),  we  find: 


Q = S7fi  | 'wg'^— '• 


This  equation  is  also  the  solution  to  our  external  problem. 

For  the  purpose  of  further  analysis,  it  is  also  necessary  to  determine 
the  non-uniform  distribution  of  heat  losses  over  the  cylinder  surface.  For 
this  goal,  as  in  solving  the  "internal"  problem,  we  introduce  the  relationship: 

' (48) 


Here  Q is  used  to  indicate  the  heat  losses  from  the  flow  to  that 
e 

half  of  the  cylinder  surface  which  is  turned  towards  the  filtration  flow. 

In  order  to  determine  the  value  Qg,  we  examine  equation  (44)  which 
in  this  case  is  rewritten  as  follows: 


Qe= 4j/ 


where  n is  found  from  (40)  when  x = 0 and  y = R : 
c o 

*)c=2tf0l'-,  (50) 

and  n retains  the  value  given  to  it  in  (45) . 

By  substituting  the  values  which  we  have  found  for  and  nc  into 
equation  (49),  we  obtain: 


Qe = 3,64  /ft  | / . 


Now,  by  solving  (47),  (48)  and  (51)  together,  we  find  the  desire 
value  of  E ' : 

E :0.  / 0'j  ---  cons*.  (5^) 

C.  Combined  Solution  to  the  "Internal"  and  "External"  Problems 


Having  examined  two  portions  of  one  problem,  we  can  now  turn  to 
obtaining  general  computation  functions. 
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In  view  of  the  fact  that  under  the  conditions  of  established 
thermal  conditions  the  quantity  of  heat  extracted  from  the  filtration  flow  is 
equal  to  the  quantity  of  heat  extracted  by  the  tube  from  a frozen  column  of 
ground,  on  the  basis  of  equations  (23)  and  (27)  we  can  write: 


— i £ll_  = s’j/  >$**4  .*? ,/.  rfi, 

— In  V ~ 


, _ i R.j  --  r-\  — e , , ^ , (54) 

j v r ~r  a/  — — 1 ; 

~t'~rr  f/  *e-ury 

o = J_  /"  [(/?,) -4)  -t-cj'’  . 

-«*•>  ‘ 3/  4,)AV;  ’ ' (55) 

After  very  simple  transformations,  equation  (53)  can  be  reduced 
to  the  following  form: 


1 }j__  -3  J 7~  = r_  ..1., 

'/Vw  V'"wg  VV  Tii  * 1 ° n ' 


In  the  terms  of  the  essence  of  the  problem,  all  of  the  values 
which  belong  to  equations  (54)  , (55)  and  (56)  except  for  Rq  and  eQ  should 

be  given.  Thus,  equation  (56)  will  contain  two  unknowns:  Rq  and  e . 

In  order  to  determine  eQ  we  can  utilize  equations  (30)  and  (52) 
from  which  it  follows  that  with  .established  thermal  conditions  E = E1  , or: 


i arc  ««(-!/ 


rT 


- ' —0,705. 

X1, 


By  transforming  this  equation  and  by  solving  it  in  terms  of  e , 
we  obtain  the  quadratic  equation: 

el-hy  l-7f5^  - Iri-e.-i  R]r~  ri.,-0.  (5g) 

From  this  equation  it  is  easy  to  determine  the  desired  e : 

o 

e{t=*l,2Rn—\/  0,44/?;,  -rr  (59) 


Together  with  (56)  this  equation  also  serves  as  the  basis  for  solving 
the  problem  about  the  boundary-equilibrium  state  of  a frozen  ground  column 
in  the  field  of  a filtration  flow. 

3.  Sample  Problem  Solution 

Assume  that  we  know: 

Temperature  on  the  surface  of  freezing  tube  T = -20°C 

Outer  radius  of  freezing  tube  r^  = 0.10  meters 

Velocity  of  filtration  flow  V = 1 m/hr 

Temperature  of  filtration  flow  T.  = +5°C 

fl  kcal 

Thermal  conductivity  factor  of  frozen  ground  = 2.0 


Thermal  conductivity  factor  of  ground  saturated  with  water  X = 1 m-hr*deg 
Volumetric  weight  of  water  Yw  = 1000  kg/m3 

1c  C3.1 

Heat  capacity  of  water  C = 1 r 

w kg* deg 

Find  R . 
o 

The  computation  equations  are  equations  (54),  (55),  (56)  and  (59)  which 
we  will  again  rewrite  here: 


4]/  1/t^g ' rfi]/V  ^ R ‘ !n  B , ’ 


' W-'V-ei]1  , 

i4n~  >.,rr  ■ y 

*3  4 

B («.i  - n)  -r  eo  , / 

r [(«:,- r-r !-<?-]-  J 

2eJ,  ' y 

<?„=!, 2 fta-]/* 0.44/?:'  4-  r\  . 


By  substituting  our  given  values  of  the  computation  values  into 
equation  (60),  we  find: 


jOK'.-o-l  [ 


1/  r?jIn 

• I 1 


] 0.37 


Examining  equation  (63),  we  find: 


e -\ -r 0,0 i 


By  substituting  value  eQ  into  equations  (61)  and  (62),  we  obtain: 


-(  l.Ji?,)-  ]/ 0.44^5  t-  o.Dl ) 2 

0,2  —0,2  j/  0 U/?5  -+■  0.01 

‘ iy?- -0.01)-  ; .2R,  - V 0.44/?; -0.01  r ' 

, 0.2l.0J  — - 0,44,0 ,i  --  0.0 1 


—0.0! ; -- ; a 0 U»;— O.ii!  j" 

2,iRj  -2.'?  ,.\f  0. 14 /?,',  - 0.0 ! 


y i 


f I R:y  - o.o i ) • r 1 1 --R  - - j/  0.  HR-  4-0.01 ' 


•>,4Ri-2R;  i/  0.44=;  -0,01 


In  view  of  the  fact  that  solving  equation  (64)  in  the  presence  of  (65) 
and  (66)  relative  to  Rq  is  impossible,  we  will  determine  the  desired  Rq  value 

by  selection.  Thus, 

when  R =0.2  meters:  | R~  !n  "-=0,225 

O ' B„ 

when  R =0.3  meters:  1 ~Rl  In  — =0,57  .«v. 


Thus,  by  comparing  the  latter  solutions  with  equation  (64),  we  see  that 


the  desired  Rq  lies: 


0,2  0?,  <0,0. 


Now  assigning  the  value  Rq  = 0.25,  we  find: 


from  which  we  may  conclude  that  the  desired  solution  will  correspond  to  the 

value  R =0.24  meters, 
o 
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velocity  of  filtration  flow; 
volumetric  weight  of  water; 
heat  capacity  of  water; 

thermal  conductivity  factor  of  frozen  ground; 
thermal  conductivity  factor  of  thawed  ground; 
radius  of  freezing  tube; 
radius  of  frozen  column  of  ground. 


This  example  shows  that  the  problem  solutions  are  associated  with  fairly 
cumbersome  calculations.  Therefore,  in  order  to  make  the  calculations  easier, 
we  have  drawn  up  the  nomogram  presented  in  Figure  3. 


In  order  to  plot  the  nomogram,  we  have  introduced  a new  variable: 


.>= 1=..  (67) 

4 y \cw  > Awg  > v 

The  temperature  oh  the  surface  of  the  freezing  tube  is  plotted  on  scale  I. 
On  scale  II  is  the  relationship  of  the  temperature  on  the  freezing  tube 
surface  to  the  temperature  of  the  filtration  flow. 

A unique  line  corresponds  to  each  value  of  the  filtration  flow  temperature 
in  the  field  of  the  first  quadrant. 

The  value  Y'Rq  In  ~ . is  plotted  on  scale  III.  As  follows  from  equation 

(60) , with  the  aid  of  a system  of  straight  lines  in  the  second  quadrant 
corresponding  to  various  values  of  C,  the  multiplication  operation  is  performed 


If  we  examine  the  system  of  equations  (60),  (61),  (62)  and  (63),  then  we 
see  that: 


r'K.l 


(68) 


Consequently,  by  assigning  various  values  of  r^  it  is  possible  to  absolutely 
determine  the  desired  radius  of  the  frozen  ground  column  R . The  solution  to 
equation  (68)  is  also  carried  out  in  the  fourth  quadrant  with  regard  to  R . 
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The  above-cited  solution  through  a numerical  example  of  the  problem  is 
repeated  in  the  nomogram,  and  the  sequential  passage  through  the  quadrants 
of  the  nomogram  is  indicated  by  a dotted  line  with  arrows  showing  the 
direction  of  the  passage.  When 


O-- = U.'JOO 

' 4 y 1000  / 1 yrl 

value  R is  equal  to  0.25  m according  to  the  nomogram,  which  corresponds 
fairly  closely  to  the  results  obtained  previously. 


